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Abstract 

We investigate various excited states of Sine-Gordon model on a strip with 
Dirichlet boundary conditions on both boundaries using a Non Linear Integral 
Equation (NLIE) approach. 
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1 Introduction 



Consider the quantum field theory of a self-interacting scalar boson (f){x,t) on a 1+1 
dimensional strip, infinite in time direction and finite in space with bulk action 
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and with the Dirichlet boundary conditions (0, t) = 0_ + -^-m^ and t) = 
0+ + m± G Z. We shall refer to this integrable theory as the Dirichlet sine- 

Gordon (DSG) model. It has several important applications ranging from condensed 
matter physics to string theory. 

The well known bulk particle spectrum of sine-Gordon, composed of solitons and 
antisolitons with topological charge 1 and —1 respectively, and, only in the attractive 
regime < (3 < -\/47r, bound states of solitons known as breathers, is complemented 
by the rich structure of boundary bound states described in [l]. The expressions 
for the bulk S-matrices can be found in [2], those for the fundamental soliton or 
antisoliton reflection matrices are given in and the general boundary bound state 
excited reflection matrices can be found in pjj. 

An important feature of DSG is the conservation of the topological charge 
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The model enjoies the discrete symmetry of the held c/) (f) + -^m and simultaneously 

(p± ^ 4>± + ("^ ^ The charge conjugation symmetry —cj) sending 
solitons into antisolitons is also guaranteed provided (f)± — > —4>± simultaneously. It 
sends Q —Q, so one can restrict attention to the study of positive Q and then act 
with this transformation to obtain states of negative Q. The periodicity allows to 
restrict the boundary parameters to the range < (t>± < ^ } 

A problem of great interest is to connect the scattering theory approach just 
mentioned to the somewhat complementary description of perturbed conformal held 
theory data. In this letter we attack this problem from the point of view of the 
nonlinear integral equation (NLIE) approach which was developed for vacuum scaling 
functions in IHEllEl and extended later to the excited states [3 [HI El CHI HI] • In this 
framework, exact scaling functions of finite size effects provide a way to investigate the 
flows from ultraviolet (UV) to infrared (IR) scales in integrable quantum field theory, 
hence building a bridge between the perturbed CFT description and the factorized 
scattering one. The NLIE for the vacuum Casimir energy was already deduced, along 
a similar line to those we shall illustrate here, some years ago in [12]. Our interest 
here is to establish a general NLIE valid for all states in the flnite size spectrum. 

^Notice that, unlike in the single boundary case of [SEIG], in the present two boundaries model 
we cannot further restrict both boundary parameters to < < ^ . 
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This allows to identify "particle" states in the scattering description with "conformal" 
states in the perturbed CFT description. 

To make this approach viable, one should start from an exact solution of a lattice 
regularization of the model. This is normally provided, in NLIE approach, by some 
2D light-cone vertex model or equivalently by some inhomogeneous ID spin chain. By 
suitable scaling limit the continuum renormalized theory can be reached. The Bethe 
ansatz equations that solve the lattice model are turned in the continuum limit into 
a NLIE that has to be considered as the basic renormalized tool to calculate the 
eigenvalues of the integral of motions of the theory on any state in finite volume. 

The homogeneous antiferromagnetic XXZ spin 1/2 model in a chain of N sites 
with lattice spacing a, coupled to parallel magnetic fields /i+ and h_ at the left and 
right boundaries respectively, has Hamiltonian 

N-l 

7^(7, h^, h_) = -JY, + alal^, + cos7«+i) + h^a^ + h^a% . (3) 

n=l 

Here a", a = x,y,z are Pauli matrices and < 7 < vr. Where convenient, we shall 
also use the equivalent parameter p = ^ — I7 whose range is < p < cxd. The 
Hamiltonian Q, as well as its inhomogeneous generalizations, can be constructed in 
a double row transfer matrix framework. For details see |13| . 

The bare continuum limit ^ 00, a ^ while Na = L remains fixed, is known 
to give, in the homogeneous case Q a massless free boson compactified on a 
circle of radius R = 1/ ^/2{7r — 7) [14i|. Among the many possible deformations of the 
Hamiltonian Q leading to sine-Gordon in the bare continuum limit, we choose the 
one introducing alternating inhomogeneities = (— 1)"6 in the sites of the chain. 
This choice has been known for sometime [13 [H] to give a correct construction of 
sine-Gordon model in the bulk in cylindrical geometry, when the appropriate scaling 
limit is chosen, with periodic or twisted boundary conditions. It is then natural to 
expect that the same construction in the presence of boundary magnetic fields h± can 
also provide an effective tool to define the renormalized DSG theory. It is worthwhile 
to recall that, as the homogeneous XXZ chain is equivalent to a 6- vertex model on 
a square lattice, this modified XXZ chain is also equivalent to a 6-vertex model, but 
- as a consequence of the introduced inhomogeneities - defined on a lattice rotated 
by 45°, i.e. on what can be thought as a Minkovski space discretized along the light- 
cone directions. This is why this construction is often referred as light cone lattice 
construction of the sine-Gordon model [11]. For 6 ^ 00, ^ 00 and a —>■ while 
L = Na is fixed, contact can be made, along lines similar to those in [T2|, with the 
lagrangean formulation of DSG model, eq.((H). The XXZ anisotropy 7 is related to 
the SG coupling /3 by = 8(7r — 7) = Details of this bare continuum limit are 
out of the scope of the present paper. 

In this preliminary letter we put our accent on the presentation of the NLIE got 
from this construction and analyze only a few bulk excited states with nonexcited 
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boundaries, in order to show consistency with expected results. The careful treatment 
of the full situation with boundary bound states would involve more delicate issues 
of analytic continuation in the boundary parameters 0± that we choose to postpone 
to a more extensive forthcoming publication for the reasons explained at the end 
of sect. 2. 

2 Bethe Ansatz and NLIE 

The Bethe Ansatz equations for the boundary XXZ chain ^ have been written by 
Alcaraz et al. and Sklyanin [in] some years ago, using an algebraic approach. It 
is straightforward to generalize them with the introduction of the alternating inho- 
mogeneities [TOij. 

Eigenvalues of the double row transfer matrix can be constructed in terms of sets 
of distinct numbers ...^'Om called roots. They are in number of M, (M < A^) and 
must satisfy the Bethe ansatz equations 

M 

(4) 

where 

sinh -(x + iun) 

sJx) = — - 

^ sinh ^(x — iz/Tr) 

and H± is defined such that h± = sin 7 cot ^{H± + 1) and we choose as fundamental 
region —p — 1 < H± < p+1. Notice that the boundary terms in the Bethe equations 
disappear when H± = 0, i.e. h± = 1 + cos 7 = he- In such case, as it was shown in 
[213121], the system becomes SLq{2) invariant. 

The antiferromagnetic vacuum turns out to be a maximal set M = ^ of real roots 
and it exists for A^ even only. In the region < 7 < vr of interest for us, and for small 
enough boundary magnetic fields, this is the true ground state of the theory. For A^ 
odd instead the states with lowest possible total spin have M = ^^f^ roots and one 
hole. However, to deal correctly with the continuum limit one has to consider both 
N even and odd sectors, like it was shown in the periodic case in ^ [101 [H]. The 
symmetry of ^ {-dj} — > {—^j}, evident from the Bethe equations, implies that only 
roots with positive real part are independent parameters characterizing a Bethe state. 
The value 'dj = is a solution of (jH) for any A^ and M. However, the corresponding 
Bethe state would vanish, so one has always to subtract this unwanted root, i.e. to 
create a hole at -d = 0. 

The domain of root distribution can be considered as a semistrip U+ of the com- 
plex 1!} plane: 

f 22 2 

U+ = <^ ^9 G C I Re^9 > , -— < Imd < — or Re-;? = , < Im^ < — > ■ 
I 27 27 27 J 
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2 

This also excludes another unwanted root at and considers only half of the imag- 
inary axis, as it should for symmetry. However, for computational purposes, it is 
often better to double this strip by mirroring all the roots 



U = <^ eC\ Re^ G M , < Im?9 < 

I 27 27 

To each root '^j associate its mirror root i^.j = — 1?^. Define the function 

^^{'&) = Tr + i\ogs^{i9) (5) 
with the oddity condition ipJ—u) = —ipJu) fixing the fundamental branch of the 

2 

logarithm. It is periodic in the imaginary direction, with period z^, and real on the 
real axis. We choose as fundamental periodicity the strip ?9 G U. Singularities of this 
function are distributed along the imaginary axis: 

Red = 0, Imd = ±7r{k{p keZ 

so that the fundamental analyticity strip is limited to llm-i?! < 7rmin(i^,p+ 1 — u). 

In terms of the function ^ the logarithm of the Bethe equations (ji)) can be 
expressed as 



[^id^j + e) + ipi/2idj - e)] + iPH+/2idj) + ^H^/2{dj) 

2M 



k=l 



Eigenvalues of the transfer matrix can be expressed in terms of its roots. In the 
following, we shall be mainly interested in the energy spectrum, for which the formula 
is 

1 ^ / d d \ 1 d 

fc=l ^ ^ 
Define for ^9 G U the so called counting function 



Z^i^) = N [(^i/2(^ + 6) + <^i/2(^ - 6)] + + ^H^/2m 

M 

- vi{^-^k)+M^)+vim (7) 

k=-M 

in terms of which the logarithm of the Bethe equations simply becomes the condition 

Z^{i3,) = 2nl, . (8) 
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The last term in (0) takes care of the fact that in the second member of ^ the product 
does not include factors with k = j. The last but one instead explicitly subtracts the 
unwanted root t^q = 0. The integers Ij play the role of quantum numbers. 

The analytic structure of Zi\i{{}) makes it convenient to classify the roots and 
related objects of the Bethe ansatz ^ as follows 

1. real roots -^k, k = 1, ...,Mfi, dk > 0: they are strictly positive real solutions of 
iD and ©; 

2. holes -dk, k = 1,...,Nh, 'dk > 0: strictly positive real solutions of (jH]) that are 
not solutions of the Bethe Ansatz ^ ; 

3. close roots ^k, k = 1, ...,Mc: complex solutions with Re-^k > and imaginary 
part in the range < llm^^^l < minvr (l,p); 

4. wide roots ^k, k = 1,...,Miy, Re-dk > 0: complex conjugate solutions with 
imaginary part 7rmin(l,j9) < llmty^l < 

For further convenience it is useful to introduce the notion of self-conjugate roots, 

2 

i.e. wide roots with Im^^^ = whose complex conjugate is the root itself, due to 
the periodicity of Bethe equations. Their number will be indicated as Msc- Also we 
shall refer to roots or holes lying on the imaginary axis as magnetic. 

The function Z^i'd) for G M is globally monotonically increasing. However, 
there may be points where locally Zj^i^d) < 0. In particular holes or roots Sj such 
that Z^{sj) < are called special holes or roots. If a special object appears, then 
there must be also two other objects (real roots or holes) with the same quantum 
number, as imposed by the global increasing monotonicity of the counting function. 
Moreover, as two roots with the same quantum number are not allowed in Bethe 
Ansatz, at maximum one object of this triple can be a root, the others are forced 
to be holes. In the following we indicate the number of specials with A^^. A special 
object should be counted both as special (i.e. in Ns) and as root or hole (i.e. in M/j 
or Nh) according to its nature. 

One may relate the numbers of various types of roots to the 3rd component of 
total spin of the system S = ^ ~ To do that, we express the asymptotics of the 
function Zj^{d) on the real axis of d where they can be compared with the counting 
of real roots and holes. As a result we get the following counting equation to be 
satisfied by any allowed root configuration 



Nh - 2Ns = 2S + Mc + 2step(p - l)Mw + step(p 



25* 



H 



p+ 1 p+ 1 



(9) 



where [a;J denotes the integer part of x and H 

Following the standard derivation as illustrated in |H1 11131221, we obtain a NLIE 
for the function Zj^i^-d). The continuum limit can be taken by sending N ^ oo and 
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a — > in such a way that L — Na remains constant. The only way to get a sensible 
NLIE satisfied by the limiting counting function ^('J?) = lim^v-^oo Zn{'&) is to admit 
that also O rescales as 

6 ~ log ■ 



ML 

where is a mass scale. We often use in the following the dimensionless scale 
parameter / = AiL. As a result of this limit procedure, one can define the NLIE on 
the continuum 

Zi-d) = 2lsmh^ + g{^\{^k}) + P {^\H+, H_) 

-2ilm / dxG{'& - x - ie) log [l - [-l)^sc(,iz{x+ie)-]^ 

where the boundary contribution is given by 

P {'d\H+, H_) = 27r / dx[F{x, H+) + F{x, H_) + G{x) + J{x)] 
Jo 

with 

G(^) ^ -e'^^ for llm^l < 7rmin(l,p) (10) 

^ ^ 27r 2 sinh fp/c cosh f A; ' ' ^ '^^ ^ ^ 

T^o, /+~dfc ,,^sinhf(p-l)A;coshf(p + l)A; it qi ^ • n ^ 

J(t?) = / — e""^ — — — H r-4r; — for \lm.i)\ < -mm(l,p) 

^ ' 27r sinh fpA; cosh p ' '2 ^ 



F{^,H) = P^e-^-^sign(ff) ""^y^t''J''P fo^ l^-^l < "l^l (H) 
^ ' J-oo 27r ^ ^ ' sinhfpA;coshfA; ' ' 2' ' ^ ^ 

The source term is given by 

^(^|{^,}) = Y,Ck[X(k){^-^k)+X(k){^ + M 



where 



X{'&) - 27r / dxG{x) (12) 
Jo 



and {"dk} is the set of position of the various objects (holes, close and wide roots, 
specials) characterizing a certain state. They are characterized by the quantization 
rule 

P 
2 



Z^i-dj) = 27rlj , Ij eZ+f^ p = Msc mod 2 
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The coefficients Ck are given by 

_ J +1 forlioles 

\ for all otfier objects 

and for any function f{'d) we define 

{fii{'d) for wide roots 

f{'d + ie) + f{'d — is) for specials 
f{d) for all other objects 

where the second determination of fi^) is defined as 

X / o\ f f{^) + ~ ZTTsignlm-i?) if p > 1 „ ,^ „, -/in 

The contribution of special objects comes from the fact that the logarithmic term 
inside the integral can go off the fundamental branch right when < 0. In this 
case, the contribution of the jump in the logarithm amounts exactly to the source 
term of a special object. For large / where the driving term dominates, monotonicity 
excludes the presence of special holes or roots. Therefore they should not be regarded 
as objects that can be added at will, but better as artefacts that appear only for 
relatively small values of /, dictated by the breakdown of analyticity of the equation 
at certain points. It is also clear from analysis done in |Hl III3 EHl 122] that the number 
Nh — ^Ng remains constant along a flow in /, and equals Nh for / sufficiently large. 

For the vacuum state containing real roots only, this equation coincides with the 
one found some years ago in [12]. Once the equation is solved for Z{{} + ie) one can 
use this result to compute the Zi^) function at any value in the analyticity strip 
|Im^9| < 7rmin(l,p), provided the function P(^9|if+, is well defined there (see 
comments below). To extend the function outside this analyticity strip one has to 
resort to the following modification of the NLIE 

Z{^) = 2/sinh,,^9 + giM{^k]) + Pii 

-2ilm / dxGiil^ - x - ie) log [l - [-1)Msc (^z{x+ie)-]^ _ 

Once Z{;d) is known, it can be used to compute the energy. It is composed of 
bulk and boundary terms whose expression can be found in ^2] and a Casimir energy 
scaling function given by 

E = M. Ck cosh(fc) -dk — M. I — sinh xQ{x) 
k J 

In the IR limit / = TVIL — > oo one can make contact between the NLIE excited 
states and the scattering theory of an underlying field theory. The integral terms in 
the NLIE and in the energy formula go as 0(e~') and can be discarded. In the case 
of a single hole placed at -^i, the NLIE becomes 
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with 



H) = 2tt / dx 



■a 







F{x H)+ l'" rffc^,fe, sinhf fcsinhf(p-l)A; 

27r sinh jpk sinh tt/c 



Interpreted, along lines of analysis very similar to those of [IDIESI, as a quantization 
rule for momentum of a particle of energy M. cosh-i?!, this yields its reflection ampli- 
tudes at both boundaries, given by 7^(^9,^±) = e'^^^'^^\ Known identities [l2| allow 
to identify such reflection matrix with the Ghoshal Zamolodchikov soliton-soliton one, 
upon suitable identification of the parameters H± and ^± in turn are connected 
to the DSG boundary parameters 0± [H] thus allowing us to relate H± and 0± as 

Notice that the periodicity 0± — > 0± + ^ reflects in the periodicity H± H± ± 
2{p + 1). By changing sign to both magnetic flelds h± — * —h± simultaneously, one 
can check that also the antisoliton reflection matrix is correctly reproduced. 

We also checked, along similar lines, that by considering suitable combinations 
of holes and non-magnetic complex roots we reproduce the correct bulk S-matrix 
as well as the expected Ghoshal-Zamolodchikov reflection matrices in some simple 
multiparticle states. 

Our analysis is limited to the case of boundary magnetic fields \ h±.\ < he = I+C0S7 
so that there is no way to accommodate any magnetic root. This corresponds to 
absence of boundary bound states, as observed in |Tn]. All bulk states then scatter 
through the unexcited Ghoshal Zamolodchikov reflection matrix. To go beyond this 
he limitation, one should modify the Fourier representation of F, i.e. consider analytic 
continuation. The hope is that such analytic continuation should naturally introduce 
source terms for magnetic roots in g{^\{^k}) in accordance to Saleur Skorik analysis of 
Bethe equations, where they found that vacuum changes by adding the maximal string 
of magnetic roots. Boundary bound states should then be obtained by removing some 
of these magnetic roots, thus creating holes on the imaginary axis. This procedure 
turns out quite cumbersome and delicate, and needs more investigation, so we choose 
to postpone the treatment of boundary bound states to a future more detailed paper 
pTT] . In the IR limit, however, by dropping the convolution term in NLIE, it is easy 
to check that one can reproduce the Mattson Dorey excited reflection matrices. 

3 UV limit and conformal theory 

In the limit / — > we make contact with the UV regime of DSG. The roots and holes 
may rescale to infinity or stay in a finite region as 

7? = ± log y or t9 = . 
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Therefore, we classify these into three types which will be denoted by indices "±, 0". 
Accordingly, the NLIE splits into left and right kink equations. Standard manipula- 
tions [HI |22| lead to the following energy formula 



E(L) ~ - f A - — 



with 
A = 



P 



P+ 1 



1 rn^ + H_ 

2 



2p 



l_2{S-S+)-2{K~K+)] + 



p+ 1 
2p 



S 



+ N (13) 



and 



N= (/+-/+-/+ -2/+ + L+ 5+ + 2ir-2ir+ -^+-2(5+)') gZ 

where = J2k ^At with A = H^C, W ... the various types of roots and holes, and 
25+ = iV+-2iV+-M+-2M+step(p-l) , L+ = sign(p - l)(M^y - M+ ) . 
The integers K and are defined through the equations 

Z+(-oo) = 2Imlog[l - (_i)A/scgi2+(-oo)j ^ ^ ^ 2t:K 

^7+(-oo|{^9+}) = 2x(oo)(5-25+)+27rir+ . 
This result should be compared with c = 1 CFT of a boson with Dirichlet confor- 



mal boundary conditions which is compactified on a circle of radius R = [2lll25] . 
The Hilbert space is composed of Heisenberg algebra representations Qm of c = 1 
CFT whose primary states |m) are created from the vacuum by vertex operators 
: e*'"^ : of ?7(1) charge 

1/*^ 1 

2 V ^^F 2 



and conformal dimensions A^ = 2^;^. All other states are created from these ones by 
applying repeatedly the creation operators 

Qra = {a^ki---a~kp\m), ki, kp G Z+} . 

For a generic state \i) G Qm the energy is given by 

7r / 1 \ ^ 



Comparing this formula with eq. (fT!?ll we see that the winding number m is related 
to the XXZ spin by m = 25". Notice that, according to this formula, the ground 
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state (m = 0) with Dirichlet boundary condition has not conformal dimension as in 
the periodic CFT: the nontrivial boundaries contribute some energy to the Casimir 
effect. 

We end this section by presenting, only in graphical form (fig. 1), a simple example 
of numerical integration of NLIE for the vacuum and few solitonic excited states. A 
detailed numerical comparison of these and other numerical data with a Truncated 
Conformal Space Approach is planned in [31j. 

E 

: 1 




1 2 3 4 5 6 7 



Figure 1: Casimir energy levels for the vacuum state (thick line) and soliton states with 
quantum number 1 = 1 (filled line), 1 = 2 (dotted line), 1 = 3 (dotted-dashed line). In this 
example p = 2.3, H+ = 2.2 and = 2.4 



4 Conclusions and perspectives 

In this letter we have obtained the NLIE governing the finite size effects for excited 
states in sine-Gordon field theory with two boundaries, each with an independent 
Dirichlet boundary condition as a continuum limit of the Bethe ansatz equations of 
the alternating inhomogeneous XXZ spin chain. Analysis of the IR and UV limits 
gives strong evidence that the underlying model is actually the DSG theory. Under- 
standing of states describing the scattering of bulk solitonic particles with Ghoshal 
Zamolodchikov unexcited boundaries are under control. 

A better understanding of the generation mechanism of boundary bound states 
for \h±\ > he as analytic continuation of the NLIE due to the singularities of the 
boundary source terms F, and of the general structure of vacuum in this case, could 
lead to a full control of the excited boundary situations too, which represents an 
achievement of crucial importance in the framework of NLIE approach. 

Also, other bulk states should be analyzed more carefully, as breather scattering 
off unexcited or excited boundaries. In this preliminary study we did not perform 
numerical analysis of NLIE. It should however be very valuable in itself and even 
better if comparable with suitable truncated conformal space approach data. 
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Finally the understanding of this Dirichlet boundary conditions case should be 
seen just as a step towards the full investigation of NLIE for general integrable bound- 
ary conditions in sine-Gordon model, whose route has been recently opened by the 
deduction of vacuum low magnetic field NLIE in starting from the Bethe ansatz 
for non diagonal boundary conditions proposed and studied in [23 EHl EHl EH] ■ 
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